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INTRODUCTION 

Many of the hardest problems in non-Abelian gauge theories, such as QCD, are of a non- 
perturbative nature. We know from the renormalization group equations that whenever 
a process involves energies of the order of (or lower than) Aqcd> the gauge coupling 
becomes strong enough to spoil the use of perturbative methods. When we move around 
the QCD phase diagram (figure 1), changing the temperature or the baryonic chemical 
potential, the situation is more complicated and, indeed, not yet fully understood. It is 
experimentally challenging to explore different regions of this diagram however such 
programs are currently underway at RHIC and shortly at the LHC. 
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FIGURE 1. The QCD phase diagram. Many corners of the diagram are conjectural. 

If we focus on the region of vanishing chemical potential, lattice simulations have 
provided hints of the existence of a critical temperature, T c ~ 170 MeV, where a (rapid) 
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crossover to a deconfined phase occurs (in the large N c limit, the crossover becomes a 
first order phase transition). This novel phase, called the quark-gluon plasma (QGP), 
is believed to have existed during the first few tens of microseconds after the Big 
Bang. For a long time this was expected to be a weakly coupled regime. However, 
this expectation was finally contrasted with data once the relativistic heavy ion collider 
(RHIC) at Brookhaven National Laboratory started to operate a few years ago. Head-on 
Au+Au collisions at center of mass energies of y/s ~ 200 GeV/nucleon, seemingly lead 
to the production of a tiny drop of QGP. After the collision the fireball thermalizes very 
quickly, expands under its own pressure and cools while expanding with a lifetime of 
order 10 fm/c. 

Indeed it appears that, against the expectations, the QGP is strongly coupled around 
the critical temperature. Strongly coupled phenomena are generally studied using lattice 
simulations. However, to date, this approach has not been suitable to apply to real- 
time dynamical processes such as those involved in this setup. An approach based 
on a string theoretical construction, the AdS/CFT correspondence, has shown to be a 
useful tool to scrutinize non-perturbative phenomena of non-Abelian gauge theories, e.g. 
those involved in QGP phenomenology. This article aims at reviewing this formalism 
in a useful manner for particle physicists wanting to be acquainted with what seems 
to be a powerful and deep tool. After reviewing the main aspects of the AdS/CFT 
correspondence, and focusing on the case of finite temperature non-Abelian gauge 
theories, we provide a list of phenomena relevant to the physics of relativistic heavy ion 
collisions. References to more in-depth explanations are provided where appropriate. 
Other useful overviews include [1,2]. 



The Large N c Limit of Gauge Theories 

The motivation for the AdS/CFT correspondence goes back to an astonishing discov- 
ery by 't Hooft [3], that large N c gauge theories should have a strongly coupled phase 
described perturbatively in terms of closed strings. To understand why one might have 
such a description we consider a simple model with only adjoint fields {e.g. gluons), 
using the so-called 'double line notation' where the propagators of the fundamental and 
antifundamental (conjugate) degrees of freedom are written using parallel lines, with 
arrows indicating the direction of flow of the color degrees of freedom. 

With this notation, Feynman diagrams become ribbon graphs, which can be classi- 
fied topologically in terms of closed Riemann surfaces, labeled L g . For this, one must 
imagine the simplest possible Riemann surface on which the double-line Feynman di- 
agram can be accurately projected. This alone does not give us any new information, 
but taking the limit such that N c — > °° while keeping the 't Hooft coupling, A = g\ M N c , 
fixed, we find that these ribbon graphs can be reordered into a series of topological 
diagrams which will be much simpler than the original Feynman graphs. In these rib- 
bon graphs, vertices (V), propagators (P) and loops (L) contribute, respectively, with 
factors of g Y ^, g\ M and N c to the amplitude of each diagram. Then, the X and N c con- 
tributions to a particular graph can be written in a factorized form as c g (X)N?, where 
X = V — P + L = 2 — 2g — Eul(L g ) is the Euler character of the given Riemman surface. 



In figure 2 we have two diagrams with the same number of vertices and propagators, one 
of which is planar, and can therefore live on a two dimensional surface without holes, 
while the other is non-planar and may be drawn on a torus faithfully. The double-line 
notation makes evident that the latter contributes a factor N% less than the former. In the 
limit that N c — > °° only the planar diagram is important. 




FIGURE 2. Two diagrams which contribute to the same order in the coupling constant, (gy^) 2 ' since 
both have 4 vertices and 6 propagators. The diagram on the left has four loops and will be proportional to 
Nc(gy M N c ) 2 , whereas the two loops on the right diagram only contribute a factor of N®(gy M N c ) 2 . 

Any amplitude for gluon interactions, si ' , can in fact be written as a sum over 
topologies 

oo 

sf=Z,c g (X)Nt 2g , (1) 

g=0 

and in the 't Hooft limit only the most trivial topology will survive. It is important to 
note that while this large N c gauge theory is clearly not the same as QCD, it has been 
shown using lattice simulations [4], that in many cases the results for large N c are almost 
the same as those for N c = 3 theories. 

We see that there are two regimes for this theory when parameterized by the coupling 
A. For A 1 the theory is perturbative in A and we can simply use Feynman diagram 
summation to calculate amplitudes. However, for A 3> 1 an increasing number of dia- 
grams contribute at each order in N c and the ribbon graphs become dense. They become 
discretized two dimensional surfaces and we can think of these as the worldsheet of 
a string theory. Therefore in the 't Hooft limit, not only do we retain only the planar 
graphs, but also these graphs are described in terms of a simple topology of a string 
world sheet. Indeed, the amplitude can be rewritten in string theory notation as 

oo 

^=l>, 2g -V g (A), g s = l/N c <l, (2) 

8=0 

where s/ g (X) is a perturbative closed string amplitude on a Riemman surface T, g , g s is 
the string coupling, and A = g s N c is a modulus of the target space. We would like to 
interpret this as a sum over topologies of the QCD string. 
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FIGURE 3. Appropriate perturbative description of a planar gauge theory as a function of X. 

From the above analysis we can see that the gauge theory admits two complimentary 
descriptions, as shown in figure 3. This gives us hints that we may be able to understand 



QCD in the strong coupling regime using string theory, but clearly we need to know 
which is the appropriate string theory to describe QCD in this limit. In order to approach 
this theory we must understand more about explicit constructions of the gauge/string 
duality. 



Historical Notes 

We sketch here some of the most important points in the history of our understanding 
of the gauge/string duality in order to comprehend with more clarity how we made such 
great leaps forward since 't Hooft's original proposal. 

In the previous section we saw how a gauge theory at strong coupling could be 
descibed by a weakly coupled string theory. However, in the early nineties, by looking 
at the case of lower dimensional matrix models, Polyakov noted that these dual string 
theories should have an extra dimension [5] on top of those describing the weakly 
coupled phase. Indeed in the case of QCD in four dimensions, he realized that this 
may account for the width of the QCD string, while the penetration depth in the fifth 
dimension would provide a natural scale. Polyakov later noted that the extra dimension 
had to be warped [6]. On its own, this was not enough to develop a concrete realisation 
of this duality in four dimensions. However as new ingredients of string theory, called 
D-branes, were discovered by Polchinski [7] the subject was dramatically altered and the 
consequences that this could have for such gauge/string dualities were quickly realised. 
This advance came when Maldacena took Polchinski's D-branes and applied Polyakov's 
ideas to come up with the first highly non-trivial four-dimensional example of a QFT 
exactly exhibiting 't Hooft's original ideas of a gauge/string duality [8]. 

THE ADS/CFT CORRESPONDENCE 

String theory is not merely a theory of strings. The other ingredients in the theory are 
now well understood (in some regimes) and include p-dimensional membranes, known 
as Dp-branes. The possible dimensionality of these objets is fixed by the string theory 
in question but here we will concentrate on the case of D3-branes, which are 3+1 
dimensional hypers urf aces. These are solutions of type IIB string theory, and admit a 
two-fold description. 

The first description is in terms of open strings. D3-branes act as hypersurfaces on 
which open string endpoints can live. One can study the low-energy description of such 
open strings, which is given by jV = 4 U (N c ) superconformal Yang-Mills (from now 
on, simply SYM) theory. 

Alternatively, from the closed string point of view D3-branes are solitonic solutions of 
type IIB supergravity. They act as sources for closed strings, whose massless spectrum 
includes the graviton and the Ramond-Ramond (RR) self-dual five-form field strength. 
The AdS/CFT correspondence emerges from the complementarity of these open and 
closed string descriptions of D3-branes. 



SYM Theory 



SYM theory in four dimensions (the dimensionality of the world-volume of the D3- 
branes) has a field content given by a superfield which is composed of one vector, A^, 

six scalars fy 1 (I = 1...6), and four fermions % x a , ( M = 1;2,3,4) which are in the 4 
and 4 representations of the SU (4) = 50(6) 7?-symmetry group. The SYM Lagrangian 
contains two free parameters g Y M and Oym, 



~^sym = Tr 



1 



2^m 



7,7 



Sx 2 



+^ x . (3) 



Not only does SYM not have any explicit scale in the classical Lagrangian, but even 
quantum mechanically since the theory is superconformal and maximally supersym- 
metric. The maximal superconformal group in four dimensions is PSU(2, 2|4) whose 
bosonic sector is 50(4,2) x 50(6). On top of this large set of global symmetries, the 
theory displays a strong/weak duality, which acts on the complexified coupling constant 
as 

9ym , 47H 1 

T ™ = — + — i 1ym^>--_ — , (4) 
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which combines with shifts in Oym to close the group 5L(2,Z). 



Multiple D3-branes in type IIB supergravity 

At low energies where the length of the string is unimportant, type IIB string theory 
is well approximated by type IIB supergravity. The supergravity solution describing a 
stack of N c coincident D3-branes is given by a metric of the form 

/ r4\-!/ 2 / t4\ V 2 

ds 2 =(l + -^\ [-dt 2 + dx 2 } + ll + -^J [dr 2 + r 2 dQ|], (5) 

plus a constant dilaton 4>, the vev of which is related to the string coupling (g s = 
e ), and N c units of Fr 5 i flux. L is the only length scale in the solution. This metric 
interpolates between a throat geometry with curvature dictated by L and an asymptotic 
ten dimensional Minkowski region. 

If we take the near horizon limit of the solution given in eq.(5), r<L, and redefine 
z = L 2 /r, we can completely decouple the Minkowski region and are left with a throat 
geometry which is given by 

[2 

ds 2 = -2 [-dt 2 + dx 2 + dz 2 ] + L 2 dQ. 2 5 , (6) 

which is the Poincare wedge of the direct product of five dimensional anti-de-Sitter 
space and a five sphere (AdSsX S 5 ). The isometry group of this space is given by 
50(4,2) x 50(6), though if we include fermions, the full supersymmetric isometry 



group is PSU(2,2\4). Note that this is exactly the same as the full global symmetry 
group of the low energy limit of the open string sector (i.e. SYM theory). All fields in the 
closed string sector come in multiplets of this supergroup, while in the open string sector 
the multiplets define BPS operators [9]. The 5L(2, Z) weak/strong coupling symmetry of 
the SYM theory manifests itself when we identify g\ M = Ang s and 6ym = 2nx, where 
X is the vacuum expectation value of the RR scalar (the axion - one of the closed string 
degrees of freedom which is classically zero for the AdSsx S 5 solution). The closed 
string sector then has precisely the same 5L(2, Z) symmetry acting on the dilaton-axion 
degrees of freedom. 

Because the AdSsx S 5 geometry is a BPS solution of the IIB supergravity equations 
of motion, one can identify the ADM mass of the D3-brane solution with N c times 
the tension of a single D3-brane (the charge), L 4 Q. 5 /4nG\o = N c \/n/K, where Gio 
is the ten dimensional Newton constant, £1$ = 7Z 3 is the volume of the 5-sphere, K = 
V / 87tGio = Sn 7 ^ 2 g s a' 2 , and a' is related to the string tension. We see that the radius L, 
of both the AdS throat and the S 5 , in string units is given in terms of the gauge theory 
parameters as: 

L 4 = g 2 YM N c a' 2 = Xa' 2 . (7) 

Therefore, in order that the stringy modes be unimportant, L 3> \ftt i , which translates 
into gauge theory language as A = gy M N c 3> 1. This exactly matches the expectations 
from 't Hooft's construction and we see that the region where the low energy super- 
gravity limit is valid is exactly where the gauge theory is non-perturbative. This is one 
the most powerful observations of the gauge/string duality and allows us to understand 
non-perturbative gauge theories in terms of low energy supergravity. It is important to 
reiterate here that by analyzing the near horizon (low energy) limit of N c coincident flat 
D3-branes, we have accumulated evidence supporting the claim that [8]: 

Planar SU (N c ) SYM theory is dual to type IIB supergravity on AdSs x S 5 
with N c — > °° units of flux. 



The AdS/CFT Correspondence 

As we have seen in the last section, when one half of the open/closed string dual 
description of D3-branes is strongly coupled the other is weakly coupled. This means 
that it is very difficult to prove (or disprove) the conjectured duality. There are indeed 
three versions of the conjecture. These are: 

• The strong version: Type IIB string theory on AdSs x S 5 Qjg s and VL 2 / a') is dual 
to SU (N c ) SYM theory (V g ¥M and V7Y C ) . 

• The mild version: Classical type IIB strings on AdSs x S 5 (g s — > and L 2 / a' fixed) 
is dual to planar SU (N c ) SYM theory (N c — > °° and A = gy M N c fixed). 

• The weak version: Classical type IIB supergravity on AdSsx S 5 (g s — > and 
L 2 /a' — > oo) is dual to planar SU(N C ) SYM theory at strong coupling (N c — > °° 
and A — > oo). 



Which of these holds is still under debate but more and more evidence is constantly 
mounting in their favor (most evidence is within the more tractable regime of the weak 
conjecture). 

In order to use the duality to perform calculations, we need a dictionary which relates 
calculations on one side with calculations on the other. The hope is then that while a 
calculation may be impossible (or very difficult) on the strongly coupled side, it may be 
very easy on the weakly coupled, dual side. To establish the dictionary between the dual 
descriptions is an ongoing and very difficult process. We will now comment on some of 
its most important entries. 



Correlation Functions 

Here we discuss the dictionary entry which relates correlation functions on either side 
of the duality. Consider a supergravity field <J» ; of mass m ; which can be shown to have 
an asymptotic solution close to the AdS boundary 

^•(z,x^)~<p ; (x^)z 4 - Ai , where A f = 2+ + (8) 

The golden entry in the holographic dictionary relates the partition function of the 
string theory to the generating functional of correlation functions of the SYM theory. 
In particular a gauge theory operator 0[ is associated with fluctuations of a (dual) 
supergravity field, 4> ; , through the following relation [9, 10] 



string 



= ^exp(y>x<p ( V^ , (9) 



where the boundary value for the supergravity field acts as the source for the field 
theory operator. The key question then is which field can source which operator. Here 
we exploit the fact that the global symmetries of the two sides of the correspondence 
match, to find a single supergravity field matching a single operator to give a product 
which is a scalar under all the global symmetries. In particular, by looking at the 
dilatation symmetry of the superconformal group we find that the the operator G\ must 
have conformal dimension A ; . Therefore for any gauge invariant operator there exists a 
corresponding string state whose mass is related to the scaling dimension of the former 
and vice-versa. 

For most applications, we are in the classical limit (g s — > 0) and it is sufficient to deal 
with the saddle-point approximation, 

iTstnng [fy] « exp (-rS5? } [<&,-]) . (10) 



Holography: the Radius/Energy duality 



Another extremely important entry in the AdS/CFT dictionary relates the radial direc- 
tion of the supergravity geometry to the energy scale in the field theory. The AdSs factor 



can be thought of as a warped (codimension one) M 1 ' 3 space, 

L 2 

ds 2 = [-dt 2 + dx 2 + dz 2 ] . (11) 

Proper times and distances in the bulk and their field theory counterparts are related 
through the warp factor by 

L z 

£SYM = - ^proper and ^SYM = ^proper- (12) 

This means that as z — > 0, we are really describing the UV of the gauge theory. Indeed 
UV divergences in the field theory are related to IR divergences in the gravity theory 
[11]. The radial direction in AdSs is described as a holographic coordinate that, from 
the gauge theory point of view, amounts to the energy scale. This is probably one of the 
deepest aspects of the AdS/CFT correspondence whose implications may be, if correct, 
as revolutionary to the understanding of spacetime dimensionality as the inclusion of 
time on an equal footing as space in special relativity. 



Wilson Loops and the qq potential 

In field theory the Wilson loop is an important tool in understanding properties such 
as confinement. However, for a strongly coupled gauge theory the calculation of the 
Wilson loop is very difficult. We may ask whether there is a way to calculate this object 
using the AdS/CFT correspondence and come up with a dual string theory description 
of the Wilson loop. In the field theory, the Wilson loop is defined as 



W[tf] = — TrP 



e i h A 



(13) 



where the integral over the gauge field is taken over some closed loop ^ and depends 
on the representation of the gauge group. In the fundamental representation, it can be 
thought of as a moving qq pair. Choosing ^ to be a rectangle, with the qq pair a distance 
L apart and propagating for a time T, the potential between the quarks, V q q(L), can be 
read off in the limit that T — > oo 5 via 

(W[<#\) =A(L)e- TY ^ L K (14) 

From the open string description of the D3-branes, a quark has a string ending on it. In 
order to study a qq pair we introduce a W boson which is associated with the breaking 
of U(N C + 1) — > U(N C ) x £7(1). From the D3-brane point of view this corresponds to 
separating one brane from the stack of N c + 1 branes. In the 't Hooft limit, the appropriate 
description of the N c branes is given by their supergravity solution and we are left with a 
single D3-brane in the AdSs x S 5 background. The qq pair is now described by a single 
string starting and ending on the remaining D3 -brane (see figure 4). Although attached 
with both ends on the brane, the rest of the string is free to explore the bulk. Given the 
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FIGURE 4. We introduce the quark pair in terms of the W bosons associated to the symmetry breaking 
U(N C + 1) — > U(N C ) x U(l), which we perform by separating a D3-brane. At strong coupling, the N c 
D3-branes are replaced by the background. 



FIGURE 5. In a theory like QCD, the string explores from the Minkowski boundary to the bulk up to a 
given radius, this effectively provides a width to the gauge theory confining string. 



separation between the ends there is a natural distance that the string probes into the 
bulk of the geometry, see figure 5. 

The expectation value of the Wilson Loop is then calculated from the partition func- 
tion of the classical string with boundary conditions specified by taken in the saddle- 
point approximation. The Nambu-Goto (NG) action gives the solution to the open string 
catenary problem [12, 13] 



In the case of SYM theory, we find the very simple result that the L dependence is 
Coulombic, as expected from the highly restrictive conformal invariance. However, the 

dependence on the 't Hooft coupling is more subtle, Y^ ong \L) = &{X) V^ eak ^(L). 
By looking at this potential it is clear that SYM is not a confining theory. As a general 
rule it is possible to determine whether the field theory dual to a given supergravity 
background is confining or not, just by minimizing the Nambu-Goto action with given 
boundary conditions. 



Up to now we have been focusing on the supergravity dual of SYM theory described 
by fields propagating in AdSsx S 5 . In this case the AdS/CFT dictionary is very well 
understood and many non-trivial tests have given agreement between the two sides of 




(W[^])~exp[-5 NG (^)]. 



(15) 



SYM vs QCD 



the correspondence [14]. Can we however think of SYM as a toy model for QCD? The 
most important differences between the two theories are: 

• QCD confines while SYM is not confining. 

• QCD has a chiral condensate while SYM has no chiral condensate. 

• QCD has a discrete spectrum while that of SYM is continuous. 

• QCD has a running coupling while SYM has a tunable coupling and is conformal. 

• QCD has quarks while SYM has adjoint matter. 

• QCD is not supersymmetric while SYM is maximally supersymmetric. 

• QCD has N c = 3 while SYM has N c -> °°. 

The answer is clearly no ! SYM and QCD are two vastly different theories and trying to 
understand QCD through SYM is not going to get us very far. However, there are ways 
that we can incrementally change the supergravity theory such that it does describe a 
better toy model of real QCD. 



Towards QCD 

It is possible to approach a QCD-like holographic dual, addressing the above issues 
in a variety of ways. Let us give a flavor of the various possibilities and some helpful 
references for the interested reader. See [15] for a recent review on this subject. 

• It is possible to consider multiple D3-branes on curved backgrounds (by replacing 
the sphere with a more generic manifold), leading to an interesting family of = 1 
superconformal field theories [16, 17]. These theories contain matter fields. 

• The theories above can be further deformed to break conformal symmetry and 
introduce confinement. This moreover leads to chiral symmetry breaking and a 
running coupling constant. These theories, in addition, display novel features like 
a cascade of Seiberg dualities along the renormalization group flow [18]. 

• Another avenue is to consider higher dimensional D-branes wrapped on certain sub- 
manifolds of the ten dimensional geometry. It is possible to obtain in this fashion 
theories that look like ,JV = 1 supersymmetric Yang-Mills theory in the IR [19, 20]. 

• Various deformations of the original geometry lead to non- supersymmetric, non- 
conformal gauge theories which exhibit confinement and chiral symmetry breaking 
[21,22,23,24,25]. 

• Fundamental matter in the quenched approximation, Nf <C N c , can be added to 
the gauge theory by introducing D7-branes [26], ignoring their effect on the back- 
ground geometry (see figure 6). Adding dynamical quarks, beyond the quenched 
approximation, is a far harder problem though much progress has also been re- 
cently made in this direction [28]. 

• Recently, a bottom-up approach has been suggested in which phenomenologically 
viable models are constructed, motivated by the AdS/CFT machinery but not within 
the full string theory framework. The link between cut-off AdS spaces and strongly 
coupled QCD may be motivated from the light-front holography approach, with 




FIGURE 6. A stack of N c D3-branes and Nf D7-branes. If Nf *C N c , the near horizon limit replaces the 
D3-branes by the background AdSs x S 5 , while the D7-branes deform due to the warped geometry. Open 
strings stretching between a D7 and a D3-brane are in the fundamental of both the color group and the 
global flavor group, while strings which have both ends on the D7 are singlets under the color group and 
adjoint under the flavor group. They source the mesons of the S YM theory [27]. 



interesting results for hadron phenomenology (see [29] and references therein). 
This line of research is known as AdS/QCD and has been quite successful in 
reproducing low energy QCD spectroscopy [30, 31]. 
• A string theory based approach similar to AdS/QCD uses so-called non-critical 
constructions in d ^ 10 dimensions [32, 33]. This makes contact with or attempts 
to motivate the AdS/QCD approach [34]. 



ADS/CFT AT FINITE TEMPERATURE 

Thermodynamics 

Introducing temperature T means adding energy without modifying other quantum 
numbers such as the brane charges. Gravity duals then involve black branes in an 
asymptotically AdS background which is simply a black hole embedded in AdSs [9]. 
The ten dimensional metric is given by 



ds 2 



r 

I? 



- 1- 



H 



dt + dx 



+ 



1- 



r 4 
JL 

A 



dr 2 + r 2 da 2 5 



(16) 



The asymptotic behavior reveals that the UV physics is not affected by the presence of 
temperature. However, the IR physics is dramatically modified. The metric above has 
a regular event horizon at r = r# with the Hawking temperature, 7#, of the black hole 
identified as that of the dual SYM theory, T = Th = rn / 7tL 2 . The Bekenstein-Hawking 
(BH) entropy density of the black hole is proportional to the area of the event horizon 
[35, 36] 

A h Tt 2 
4hG " BH = T 



-Af c 2 r 3 . 



(17) 



Notice that the entropy, besides being proportional to N 2 (a signal of deconfinement), is 
3/4 of the Stefan-Boltzmann value for a free gas of quarks and gluons. This factor is not 



a mistake of the AdS/CFT calculation, rather it is a prediction for a finite temperature 
strongly coupled SYM theory. It appears that 

*sym = ^(A)s S b, ^(0) = 1, &H = \- (18) 

Indeed, perturbative computations [37, 38] and string theory corrections (at large X) 
[39], are consistent with a monotonic function &{X). 

Lattice calculations of the energy density of QCD displays a plateau with e ~ 4/5 £sb 
[40]. Notice that 4/5 is closer to 3/4 than to 1. The result for SYM loosely suggests that 
the value of the entropy density in lattice simulations might be interpreted as a fingerprint 
of a strongly coupled quark-gluon plasma. The question now is whether or not SYM at 
finite temperature is closer to real QCD than its zero temperature counterpart. 



SYM vs QCD at Finite Temperature 

The AdS/CFT dictionary is of great help in understanding SYM theory at finite 
temperature. Just as we did for zero temperature, we can make a comparison between 
SYM and QCD in the heated phase: 

• None of them is confining, for T > T c . 

• None of them is conformal. 

• QCD is strongly coupled at T c <T < 2T C , SYM has a tunable coupling. 

• Both theories display Debye screening. 

• Both theories are non supersymmetric. 

The finite number of colors and the inclusion of matter in the fundamental representation 
are clearly not altered by the finite temperature set up and these should be separately 
addressed as before. 

We can see clearly that while the theories are still not identical, there are much 
stronger similarities than in the zero temperature case and we may hope to gain some 
more insight into real world finite temperature QCD using the AdS/CFT correspondence 
than we could do at zero temperature. 

ADS/CFT AND THE QUARK-GLUON PLASMA 

As discussed in the introduction, the experiment at RHIC is an excellent testing ground 
which we can use to compare real-world QCD with the results from the AdS/CFT 
correspondence at finite temperature. There are many properties of the strongly coupled 
plasma which can be studied but we will concentrate on just a few here, leaving a short 
discussion of several others for the end. 



Elliptic Flow 



At heavy ion colliders, in off-center collisions, (impact parameter b ^ 0), the heated 
overlap region is almond shaped. Either free streaming washes out the 'almond', or else 
collective interactions result in an anisotropy of detected particles with respect to the 
interaction plane. From such collisions (the distribution of which is well understood) 
we can ask what the number of particles detected at a given angle will be. This is 
parameterized by 

dN 

— « l + v 2 (A,£, ft )cos20, (19) 

where the parameter v 2 , the second moment, which depends on the atomic number of 
the colliding ions, A, the impact parameter, and the transverse momenta of the particles, 
p t , tells us about the elliptic flow within the medium after the collision. Hydrodynamics 
seemingly provides a good description of the QGP [41]. The experimental results for v 2 
can be compared with models which are parameterized as a function of r\/s, the shear 
viscosity to entropy ratio. The results of such a comparison are sketched in figure 7. It is 
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FIGURE 7. Behavior of V2 as a function of p t (A and b are kept fixed) for varying, small values of 77 /s. 
Hydrodynamics is destroyed for high values of b, i.e. highly off-center collisions. 

found that the best fit to the data is obtained for rj/s = 0.15 h/ kg and is certainly below 
5/(4n)h/kB [42]. The perturbative value of r\/s goes to °o for X — > [43], which is 
manifestly in conflict with RHIC data, therefore it would be useful to scrutinize its value 
in the strongly coupled regime. 

The calculation of the shear viscosity comes from studying the finite temperature 
stress energy tensor which, in the long wavelength limit, is of the form 



Tij = p 8, 



v 



r\ ^diUj + djUi-^8ijd k Ukj , 



(20) 



where p is the pressure and w, is the local fluid velocity. Using the Kubo formula it can 
be shown that [44] 



tj = lim -L / dtd'xe^HT^t^), ^(0,0)}) 

(0^0 2(1),' 



(21) 



an expression that can readily be put through the AdS/CFT machinery. Recall that the 
golden entry of the AdS/CFT correspondence, eq.(9), allows us to compute just such 
a two point function. In fact, the calculation amounts to finding the absorption cross 



section, a a b s , of a low energy graviton h xy (the field which couples to the stress energy 
tensor on the boundary) by the black hole, rj = a a b s ((U — > 0)/\6nG [45, 46]. Under 
very general conditions, C7 a bs(<y — ► 0) = An [47, 48]. This, in turn, is proportional to the 
entropy density. Then [44] 

n = *-N^ =► ^- = ~. (22) 
' 8 c s 47tk B v 

Based on universal properties of black hole horizons, this result appears to be extremely 
robust. All gauge/gravity duals fulfill this relation, regardless of their field content and 
number of supersymmetries. 

It is tempting and natural, then, to conjecture that it may be valid for QCD (up to 
finite A and N c corrections). Going beyond supergravity by including stringy corrections 
both to r\ and s, is consistent with a monotonic interpolation between the perturbative 
and non-perturbative results [39, 49]. This is one of the most important, and perhaps 
surprising results of the attempt to formulate a finite temperature holographic dual of 
QCD. It is, indeed, strikingly close to the experimental value and is one of the main 
successes of this approach to the physics of the QGP. 



Jet Suppression 

Another interesting signature of the QGP that one may study from the AdS/CFT 
perspective is that of jet suppression. If we consider the ratio, R, of the number of jets 
in heavy ion collisions (e.g. Au+Au) to that seen in p+p collisions (scaled to account for 
the number of participating nucleons) we see that any departure from R = 1 indicates 
that partons produced in hard scattering events are slowed by the hot medium which is 
only present in the Au+Au collisions. Indeed such suppression caused by this medium 
induced slowdown is detected at RHIC. 

The main channel for this energy loss is medium induced gluon radiation (figure 8). As 




FIGURE 8. Medium induced gluon radiation. A quark with energy E emits a gluon with a fraction of 
momentum x and transverse momentum p t . In the eikonal approximation, E ^> xE 3> Pt- 

it travels through the medium, in the eikonal approximation the particle trajectories can 
be written as Wilson lines in light-cone coordinates [50] . In the multiple soft scattering 
approximation, the transverse position of the gluon follows a Brownian motion, and 
the average transverse momentum after traveling a distance L is characterized by the 



transport coefficient, q ~ (pj)/L, called the jet quenching parameter (see [51, 52] for a 
thorough explanation). 

Calculating the cross-section for gluon emission [53] one finds that up to logarithmic 
corrections in the small distance limit 



The expectation value of the two point function for the adjoint Wilson loop and the 
corresponding medium properties are therefore all encoded in this single parameter, q. 
This expression can be extrapolated to the strong coupling limit and considered as a 
non-perturbative definition of q. Written in this limit we have the following definition 
in terms of the expectation value of the Wilson loop defined on a closed path ^ with a 
large light like side L and a much smaller space-like separation L <^ L : 



As we have seen earlier, we are able to efficiently calculate the expectation value of 
a Wilson loop using the AdS/CFT machinery. Note that in the large N c limit we are 
able to approximate accurately the adjoint Wilson loop by its fundamental counterpart: 
(W A (tf)) = (W F (tf)) 2 + 0(\/N c ) and therefore we can use the previous definition for 
the fundamental Wilson loop in terms of the NG action of the fundamental string. 

Performing the Wilson loop calculation using the AdS/CFT prescription, this time in 
the finite temperature background, we use eq.(24) to obtain <?sym 06 T 3 \/X [54]. In con- 
trast to the weak coupling expression, at strong coupling this does not depend explicitly 
on the number of degrees of freedom. We may see if this expression agrees with exper- 
imental data by considering a representative choice of values for the free parameters: 
A = 6k (i.e. a s = 1/2) and T = 300 MeV. In this case, q SYM = 4 .48 GeV 2 /fm, which 
is reassuringly close to the experimental value q exp « (10 ±5) GeV 2 /fm. It is possible 
to perform the same calculation in different geometries, corresponding to different finite 
temperature field theories. In contrast to the ratio rj /s, there is no universal value for q. 
Finite 't Hooft coupling corrections tend to reduce q, in the direction of the perturbative 
results, as one would expect for a monotonic function of A [55]. 

As explained previously, we do not have a string dual of QCD. However, we can 
try to understand how different quantities depend on supersymmetry, dimensionality, 
field content, etc. Some attempts to extrapolate results from SYM towards QCD have 
been explored recently in theories without fundamental degrees of freedom. In order to 
compare these theories with QCD one should chose an unambiguous quantity to fix in 
both theories in order to get an accurate comparison. It was suggested that the energy 
density was an appropriate quantity [56]. Choosing this comparison we are led to the 
following relationship between the SYM and QCD temperatures Tsym = 3~ 1//4 Tqcd 
which would mean that <?qcd < <7sym [57] in contrast to the above first approximation. 

It has however been alternatively conjectured that, since the QGP of QCD is approxi- 
mately conformal at T m 2T c [58], 




(23) 



{W A (tf)} ~exp 



1 



qL-L 2 . 



(24) 




(25) 



which would agree with the above tendency. Clearly the comparison is an area of 
controversy and the case is still open for discussion. It should also be noted that the 
addition of flavor in non-critical set-ups tends to challenge the latter result, suggesting 
that <7qcd may be higher than <?sym [59]. 

FURTHER FEATURES OF THE QGP 

There are many subtleties and controversial points that we do not have space to explore 
in detail here. However we can give some basic ideas of some of the most interesting 
recent investigations into using the AdS/CFT correspondence to understand several fea- 
tures of the strongly coupled QGP. We provide references to some appropriate literature. 



Meson melting 

As mentioned, it is possible to introduce fundamental matter into the AdS/CFT corre- 
spondence. Doing so allows us to study meson phenomenology and this is particularly 
interesting at finite temperature. Flavor D7-branes can be embedded into the black D3- 
brane background in two distinct ways, which give very different effects for the associ- 
ated mesons [24, 60]. The D7-brane can either end on the black hole horizon or it can 
miss the horizon completely. The solutions which do not hit the horizon (corresponding 
to large values of m q /T) display a discrete spectrum of mesons with a mass gap related 
to the value of m q . These mesons are completely stable and their spectral function is an 
infinite sum of delta functions. The solutions that touch the horizon (corresponding to 
small values of m q /T) on the other hand display a continuous spectral function and the 
associated mesons have a finite lifetime [61]. By studying this spectrum of quasinormal 
modes one can calculate the lifetime of the mesons once they enter into the QGP (or, 
more accurately the SYM plasma as there are no dynamical quarks). In particular it can 
be seen that while mesons made of light quarks melt quickly into the plasma, heavy 
quark mesons may remain relatively stable. Indeed it has been noted that J mesons 
do not appear to melt in the RHIC plasma as do the light-quark mesons. 



Photoproduction 

Photoproduction is an excellent probe of the strongly coupled quark gluon plasma. 
Once photons are produced, they travel virtually unimpeded through the optically thin 
medium and can be detected to give an idea of the interior of the hot, dense plasma. 
Although using the AdS/CFT correspondence we have no gravity dual of a weakly 
coupled U(l) theory, using the optical theorem we can write the photoproduction rate 
in terms of a two point function of vector currents, qy^q, which can be calculated by 
looking at a U(l) vector field on a probe flavor D7-brane [62, 63, 64]. Dynamical photon 
corrections will come at next to leading order to this calculation and therefore we can 
get a good handle on the true photoproduction rate without a dynamical U(l). The two 
point function can be written directly in terms of the spectral function calculated from 



excitations on the D7-brane. The effect of finite baryon density on the photoproduction 
rate can be calculated and from the spectral function we can study the conductivity and 
diffusion in the plasma. All of these calculations give more results which can in theory 
be tested against experiment. One of the major difficulties is that such experiments have 
a great deal of noise although there are some predicted signals which may be easily seen 
as sharp peaks in the spectral function [65]. 



Bulk viscosity 

The bulk viscosity is related to the energy loss of a fluid due to its expansion. For a 
conformally invariant system, £ = 0. This is the result, for instance, obtained using the 
AdS/CFT correspondence for SYM. However, as discussed throughout this lecture, there 
are other gauge/gravity systems where conformal symmetry is not apparent, and they 
yield a non-vanishing bulk viscosity. An interesting relation among different transport 
coefficients, 




where c s is the speed of sound, was obtained within this framework [66]. There is a 
family of gauge theories coming from D-brane configurations that saturate this bound 
[67], a saturation that seems to be a consequence of their kinematical structure [68]. 



Drag force and diffusion wakes 

A high energy particle moving through a hot plasma will lose energy to the surround- 
ing medium, leading to an effective viscous drag on its motion. Within the framework of 
the AdS/CFT correspondence this is represented by a string hanging from the boundary 
of AdS into the bulk. The point of attachment carries a fundamental charge under the 
gauge group SU(N C ), is infinitely massive (when the string endpoint is placed on the 
boundary) and moves with constant velocity relative to the rest frame of the plasma. 
The drag force that the trailing string inside the plasma exerts on the external quark is 
[69,70,71] 

dp n \lsYM N ^ v 

In contrast to quarks, spinning mesons, i.e. color singlet states, do not experience any 
drag effect when propagating through the plasma [72, 73]. This means that no force is 
necessary to keep them moving with fixed velocity. The same applies to baryons [74]. 

Energy transferred to the plasma from the moving quark will cause the energy density 
of the plasma, in the vicinity of the quark, to deviate from its equilibrium value. That is, 
the moving quark will create an energy density "wake" which moves with it through the 
plasma. Evaluating the energy density perturbation we observe a net surplus of energy in 
front of the quark and a net deficit behind. This may naturally be interpreted as plasma 
being pushed and displaced by the quark, just like the water displacement produced 



by a moving boat. Another interesting feature revealed by the gravity calculation is the 
formation of a conical energy wake, or sonic boom, for velocities greater than the speed 
of sound [75, 76]. Contrary to quarks where the diffusion wake is strong, when mesons 
move through the plasma the wake is absent whilst a shock wave is created [77]. 

Early times after the collision 

One may wonder whether a full description of an ultra-relativistic heavy ion collision 
may be achieved within AdS/CFT. This line of thought was advocated in [78]. The gen- 
eral approach consists of the study of black hole formation in AdS due to the collision 
of shock waves. These shock waves completely stop shortly after the collision, which 
has been interpreted as the holographic counterpart of the nuclear stopping due to strong 
coupling effects [79]. The process of horizon formation in AdS due to time-dependent 
perturbations on the boundary, oriented to the study of far-from-equilibrium isotropiza- 
tion in the plasma, was recently explored by numerical methods [80]. Moreover, by 
means of so-called holographic renormalization, it is possible to extract information on 
the real-time dynamics of the energy-momentum tensor of the gauge theory just after the 
impact (e.g. an estimate of the thermalization time as a function of the energy density) 
from the collision of shock waves [81]. 

Bjorken hydrodynamics 

The drop of quark-gluon plasma in the experimental setup undergoes a rapid expan- 
sion. From the field theory point of view, the condition of boost invariance in the cen- 
tral rapidity region, together with some other mild assumptions, leads to the specific 
pattern of Bjorken flow [41]. A holographic dual of this phenomenon was unveiled in 
[82], where the late-time dynamics reproduces the T(z) ~ t~ 1//3 cooling law. Moreover, 
studying the conditions that are necessary to avoid the appearance of singularities, the 
famous ratio r\/s=\f (An) was obtained in [83]. This analysis can be extended to higher 
orders in a long wavelength expansion, and this was shown to provide a powerful tool to 
extract all kinds of transport coefficients [84, 85]. 

There is clearly more to be discussed in this new and exciting arena, however here we 
simply refer the interested reader to the quickly growing library of research papers. 

ACKNOWLEDGMENTS 

JDE wishes to thank the organizers of the XIII Mexican School of Particles and Fields 
and specially Alberto Giiijosa and Elena Caceres for inviting him to deliver his lectures 
in a wonderful and stimulating ambience. This work is supported in part by MICINN 
and FEDER (grant FPA2008-01838), by Xunta de Galicia (Consellena de Educacion 
and grant PGIDIT06PXIB206185PR), and by the Spanish Consolider-Ingenio 2010 
Programme CPAN (CSD2007-00042). JDE is a Ramon y Cajal Research Fellow. JPS is 



supported by a Juan de la Cierva fellowship. The Centro de Estudios Cientificos (CECS) 
is funded by the Chilean Government through the Millennium Science Initiative and the 
Centers of Excellence Base Financing Program of Conicyt. CECS is also supported by 
a group of private companies which at present includes Antofagasta Minerals, Arauco, 
Empresas CMPC, Indura, Naviera Ultragas and Telefonica del Sur. 

REFERENCES 

1. R. C. Myers, and S. E. Vazquez, Class. Quant. Grav. 25, 114008 (2008), 0804.2423. 

2. S. S. Gubser, and A. Karch (2009), 0901.0935. 

3. G. 't Hooft, Nucl. Phys. B72, 461 (1974). 

4. M. Teper (2004), hep-th/0412005. 

5. A. M. Polyakov, Phys. Lett. B103, 207-210 (1981). 

6. A. M. Polyakov, Nucl. Phys. Proc. Suppl. 68, 1-8 (1998), hep-th/ 97 1 10 02. 

7. J. Polchinski, Phys. Rev. Lett. 75, 4724^1727 (1995), hep-th/9510017. 

8. J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231-252(1998), hep-th/ 97 112 00. 

9. E. Witten, Adv. Theor. Math. Phys. 2, 253-291 (1998), hep-th/9802150. 

10. S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, Phys. Lett. B428, 105-114 (1998), hep-th/ 
9802109. 

11. L. Susskind, andE. Witten (1998), hep-th/ 9805114. 

12. S.-J. Rey, and J.-T. Yee, Eur. Phys. J. C22, 379-394 (2001), hep-th/9803001. 

13. J. M. Maldacena, Phys. Rev. Lett. 80, 4859^862 (1998), hep-th/ 9803002. 

14. I. R. Klebanov, AIP Conf. Proc. 1031, 3 (2008). 

15. J. D. Edelstein, and R. Portugues, Fortsch. Phys. 54, 525-579 (2006), hep-th/0602021. 

16. S. S. Gubser, Phys. Rev. D59, 025006 (1999), hep-th/9807164. 

17. I. R. Klebanov, andE. Witten, Nucl. Phys. B536, 199-218 (1998), hep-th/ 98 07 80. 

18. I. R. Klebanov, and M. J. Strassler, JHEP 08, 052 (2000), hep- th/ 07191. 

19. J. M. Maldacena, and C. Nunez, Phys. Rev. Lett. 86, 588-591 (2001), hep- th/ 00 8 001. 

20. J. D. Edelstein, and C. Nunez, JHEP 04, 028 (2001), hep-th/0103167. 

21. N. R. Constable, and R. C. Myers, JHEP 11, 020 (1999), hep-th/ 9 9 05 81. 

22. J. Polchinski, and M. J. Strassler, The string dual of a confining four-dimensional gauge theory 
(2000), hep-th/ 0003 13 6. 

23. J. Polchinski, and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002), hep-th/010 9174. 

24. J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik, and I. Kirsch, Phys. Rev. D69, 066007 (2004), 
hep-th/0306018. 

25. T. Sakai, and S. Sugimoto, Prog. Theor. Phys. 113, 843-882 (2005), hep-th/0412141. 

26. A. Karch, and E. Katz, JHEP 06, 043 (2002), hep-th/0205236. 

27. M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, JHEP 07, 049 (2003), hep-th/ 
0304032. 

28. R. Casero, C. Nunez, and A. Paredes, Phys. Rev. D73, 086005 (2006), hep-th/ 6 02 02 7. 

29. S. J. Brodsky, and G. F. de Teramond (2008), 0812.3192. 

30. J. Erlich, E. Katz, D. T. Son, and M. A. Stephanov, Phys. Rev. Lett. 95, 261602 (2005), hep-ph/ 
0501128. 

31. L. Da Rold, and A. Pomarol, Nucl. Phys. B721, 79-97 (2005), hep-ph/0501218. 

32. A. M. Polyakov, Int. J. Mod. Phys. A14, 645-658 (1999), hep-th/9809057. 

33. I. R. Klebanov, and J. M. Maldacena, Int. J. Mod. Phys. A19, 5003-5016 (2004), hep-th/ 
0409133. 

34. U. Gursoy, E. Kiritsis, L. Mazzanti, andE Nitti, Phys. Rev. Lett. 101, 181601 (2008), 0804 . 0899. 

35. J. D. Bekenstein, Phys. Rev. D7, 2333-2346 (1973). 

36. S. W. Hawking, Nature 248, 30-31 (1974). 

37. A. Fotopoulos, and T. R. Taylor, Phys. Rev. D59, 061701 (1999), hep-th/9811224. 

38. M. A. Vazquez-Mozo, Phys. Rev. D60, 106010 (1999), hep-th/ 9 905030. 

39. S. S. Gubser, I. R. Klebanov, and A. A. Tseytlin, Nucl. Phys. B534, 202-222 (1998), hep-th/ 
9805156. 



40. P. Braun-Munzinger, and J. Stachel, Nature 448, 302-309 (2007). 

41. J. D. Bjorken, Phys. Rev. D27, 140-151 (1983). 

42. H. Song, and U. W. Heinz, Phys. Rev. C78, 024902 (2008), 0805.1756. 

43. S. C. Huot, S. Jeon, and G. D. Moore, Phys. Rev. Lett. 98, 172303 (2007), hep-ph/0 60 80 62. 

44. G. Policastro, D. T. Son, and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001), hep-th/ 
0104066. 

45. S. S. Gubser, and I. R. Klebanov, Phys. Lett. B413, 41^8 (1997), hep-th/9708005. 

46. S. S. Gubser, I. R. Klebanov, and A. A. Tseytlin, Nucl. Phys. B499, 217-240 (1997), hep-th/ 
9703040. 

47. S. R. Das, G. W. Gibbons, and S. D. Mathur, Phys. Rev. Lett. 78, 417-419 (1997), hep-th/ 
9609052. 

48. R. Emparan, Phys. Rev. D56, 3591-3599 (1997), hep-th/9704204. 

49. A. Buchel, J. T. Liu, and A. O. Starinets, Nucl. Phys. B707, 56-68 (2005), hep-th/04 62 64. 

50. R. Baier, Y. L. Dokshitzer, A. H. Mueller, and D. Schiff, Nucl. Phys. B531, 403-425 (1998), 

hep-ph/9804212. 

51. J. Casalderrey-Solana, andC. A. Salgado, Acta Phys. Polon. B38, 3731-3794 (2007), 0712.3443. 

52. J. D. Edelstein, and C. A. Salgado, AIP Conf. Proc. 1031, 207-220 (2008), 0805.4515. 

53. U. A. Wiedemann, Nucl. Phys. B588, 303-344 (2000), hep-ph/0005129. 

54. H. Liu, K. Rajagopal, and U. A. Wiedemann, Phys. Rev. Lett. 97, 182301 (2006), hep-ph/ 
0605178. 

55. N. Armesto, J. D. Edelstein, and J. Mas, JHEP 09, 039 (2006), hep-ph/ 60 62 4 5. 

56. S. S. Gubser, Phys. Rev. D76, 126003 (2007), hep-th/0611272. 

57. S. S. Gubser, Nucl. Phys. B790, 175-199 (2008), hep-th/0612143. 

58. H. Liu, K. Rajagopal, and U. A. Wiedemann, JHEP 03, 066 (2007), hep-ph/ 6 12 16 8. 

59. G. Bertoldi, F. Bigazzi, A. L. Cotrone, and J. D. Edelstein, Phys. Rev. D76, 065007 (2007), hep-th/ 
0702225. 

60. I. Kirsch, Fortsch. Phys. 52, 727-826 (2004), hep-th/0406274. 

61. C. Hoyos-Badajoz, K. Landsteiner, and S. Montero, JHEP 04, 031 (2007), hep-th/ 612 169. 

62. S. Caron-Huot, P. Kovtun, G. D. Moore, A. Starinets, and L. G. Yaffe, JHEP 12, 015 (2006), 
hep-th/0607237. 

63. D. Mateos, and L. Patino, JHEP 11, 025 (2007), 0709.2168. 

64. J. Mas, J. P. Shock, J. Tarrio, and D. Zoakos, JHEP 09, 009 (2008), 0805.2601. 

65. J. Casalderrey-Solana, and D. Mateos (2008), 0806.4172. 

66. A. Buchel, Phys. Lett. B663, 286-289 (2008), 0708.3459. 

67. J. Mas, and J. Tarrio, JHEP 05, 036 (2007), hep- th / 7 030 93. 

68. I. Kanitscheider, and K. Skenderis (2009), 0901.1487. 

69. C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz, and L. G. Yaffe, JHEP 07, 013 (2006), hep-th/ 
0605158. 

70. S. S. Gubser, Phys. Rev. D74, 126005 (2006), hep-th/0605182. 

71. E. Caceres, and A. Guijosa, JHEP 11, 077 (2006), hep-th/0605235. 

72. K. Peeters, J. Sonnenschein, andM. Zamaklar, Phys. Rev. D74, 106008 (2006), hep-th/ 60 61 95. 

73. M. Chernicoff, J. A. Garcia, and A. Guijosa, JHEP 09, 068 (2006), hep-th/0 607 8 9. 

74. M. Chernicoff, and A. Guijosa, JHEP 02, 084 (2007), hep-th/ 61 11 55. 

75. S. S. Gubser, S. S. Pufu, and A. Yarom, Phys. Rev. Lett. 100, 012301 (2008), 0706.4307. 

76. P. M. Chesler, and L. G. Yaffe, Phys. Rev. Lett. 99, 152001 (2007), 0706.0368. 

77. S. S. Gubser, S. S. Pufu, and A. Yarom, JHEP 07, 108 (2008), 0711.1415. 

78. E. Shuryak, S.-J. Sin, and I. Zahed, J. Korean Phys. Soc. 50, 384-397 (2007), hep- th/ 511199. 

79. J. L. Albacete, Y. V. Kovchegov, and A. Taliotis, JHEP 07, 100 (2008), 0805.2927. 

80. P. M. Chesler, and L. G. Yaffe (2008), 0812.2053. 

81. D. Grumiller, and P. Romatschke, JHEP 08, 027 (2008), 0803 . 322 6. 

82. R. A. Janik, and R. B. Peschanski, Phys. Rev. D73, 045013 (2006), hep-th/0512162. 

83. R. A. Janik, and R. B. Peschanski, Phys. Rev. D74, 046007 (2006), hep-th/ 60 61 4 9. 

84. S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and M. Rangamani, JHEP 02, 045 (2008), 0712 . 
2456. 

85. R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, and M. A. Stephanov, JHEP 04, 100 (2008), 
0712.2451. 



